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INTRODUCTION 


Research is currently underway at Ames Research Center on developing trajectory control 
logic for a class of advanced aircraft known as vertical or short takeoff and landing aircraft 
(V/STOL) (ref. 1). From that work, simulations of such aircraft oriented to autopilot research are 
introduced as an example of a large-scale system. 

The scope of this research is to proceed from basic theory to flight test on at least one research 
STOL aircraft now at Ames. Work is currently concerned with the augmentor wing jet STOL 
research aircraft (AWJSRA); this will be the source of some of the details of aircraft models 
discussed later. Since generalized digital flight-control hardware system (called STOLAND, ref. 2) 
has been designed for Ames and installed in the AWJSRA, the results of our research can be entered 
as control logic software in the STOLAND system and flight tested. The research is focused on 
controlling the vehicle during approach and landing in anticipation of an operational environment 
of high-density terminal area traffic where tight control over flight paths, both in space and time, is 
required. 

Scale, in this paper, is sensed from the point of view of the autopilot research, that is, the size 
of the simulation is defined in relation to its use in the research. One measure of size significant to 
us was the ease with which the simulation could be used as a research tool by one or two persons. 
Another measure was its complexity as a system to be controlled. A generalized autopilot is a device 
for controlling the aircraft flight path throughout its flight regime and, as such, necessarily solves 
the aircraft equations of motion. If that solution is very difficult, then it reflects the complexity of 
the aircraft itself. 

In our research, problems of both size and complexity were encountered and these derive 
mostly from the nature of the V/STOL aircraft itself. For the class of V/STOL aircraft, several new 
design elements contribute to problems of size and complexity; they are designed to use both 
engine power and major changes in structural configuration to augment the aerodynamic forces to 
achieve low landing speeds, and this results in nonlinear engine and configuration-dependence of 
aerodynamic modeling data. Often there is also a novel type of engine control not found on 
conventional aircraft; for example, the engines of the AWJSRA exhaust through a controllable 
nozzle that can be rotated more than 90°. This control is simple to simulate but adds a great deal of 
complexity to the flight-control problem. 

The appropriate level of detail in simulating the flight dynamics is the rigid body motion of the 
aircraft which, by itself, is not a large number of degrees of freedom. 
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Problems of size as a working tool were attacked by use of a hierarchy of simulations that 
represent the aircraft with increasing completeness of detail in parallel with a corresponding hier- 
archical development of the autopilot, and also through various partitionings of the system equa- 
tions. The control problem was also attacked through partitioning to divide it into problems half as 
large although still difficult to solve. 
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SYMBOLS 

matrix mapping angular velocity into Euler attitude rates 
matrices of wind gust state equations 
matrices of approximate attitude state equation 

Ust atmospheric parameters: density, temperature, pressure, and speed of 
sound 

(^,0,'P), Euler attitude angles 
wing span 
wing cord 

dimensionless coefficients of drag, Y, and lift force 
dimensionless coefficient of equivalent thrust 
dimensionless coefficients of roll, pitch, and yaw moments 
dimensionless gradients of aerodynamic force and moment 
transformations for rotations about the axes, respectively 
dimensionless force and moment vectors 
flap deflection angle 

forces — total, aerodynamic, engine, respectively, in body axes 

aircraft inertia matrix in body axes 

inertial reference frame (runway) 

body axes reference frame 

wind axes reference frame 
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location of engine in body axes 
aircraft mass 
engine fuel flow 

moments in body axes — total, aerodynamic, and engine moments, 
respectively 

noise 

dynamic pressure 

position, vector, inertial coordinates 
wing area 
engine thrust 

equivalent thrust of compressed engine air used to augment aerodynamics 

control variable 

moment controls (5^, 6g, Sp) 

force controls (F, 6th, i') 

airspeed 

velocity vector with respect to air mass, body axes components 

inertial velocity, inertial axes components 

total, steady, gust winds, respectively; inertial axes components 

translational motion variables (Rp Vj, Vj) 

translational variables for reference trajectory 

rotational motion variables (^>, 

angle of attack, angle between body axis /g and velocity vector 
side-slip angle, angle between velocity vector and aircraft plane of symmetry 
aileron deflection angle 
elevator deflection angle 
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Subscripts 

B 

C 

EX 

I 
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Superscript 

(') 


angle variable, one of (3, 5^, 6^, 5^ 
rudder deflection angle 
throttle setting 

Euler attitude angles; roll, pitch, and heading angles, respectively 
angular rate variable, one of |a, /3, 
nondimensionalizing factor for /ay 
time constant 

angular velocity, body axes components 


quantity is a vector given in body axes 
commanded value of the quantity 
extremes of the quantity (maximum or minimum) 
quantity is a vector given in inertial axes components 
quantity is a vector given in wind axes components 

time derivative of the quantity 

SIMULATION USES IN THE RESEARCH 


The digital simulation of the system serves as a tool for both computing and understanding 
during the process of solving the control problem. First, it contains a reference definition of the 
aircraft in the form of parameters and model data and functional relationships, as well as algo-, 
rithmic specifications of the dynamic behavior of the aircraft algorithmically. All aspects of the 
aircraft influence the autopilot development so that it is advantageous for the researcher to become 
intimately familiar with the simulation. To promote these uses, the simulations should be easy to 
access, use, understand, and modify, and its structure should reflect the autopilot researcher’s 
perception of aircraft flight at as many levels as possible. Secondly, aircraft motion is sufficiently 
complex that the autopilot design begins with ideas based on a simplified view of the aircraft 
followed by an iterative use of the simulation in which ideas are tested and unforeseen effects 
uncovered • and then fed back into the design. Many effects understood only qualitatively are 
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quantitatively tested. During this process, the autopilot changes frequently while the remainder of 
the simulation program remains the same so that it will be useful to isolate the autopilot as a plug-in 
subsystem. Thirdly, the autopilot research from theory to flight test proceeds in stages of ever 
increasing integration with the complete flight-control problem and a hierarchy of simulations that 
correspondingly simulate the aircraft in increasing detail is necessary. In the first stage, since we are 
concerned with trajectory control, only a simplified simulation of the rotational dynamics is neces- 
sary. Later, attitude control and navigation is added, which requires simulation of both the transla- 
tional and rotational degrees of freedom. A final stage of integration with flight hardware, displays, 
safety, and pilot opinion is required, and a large general purpose, man-in-the-loop laboratory simula- 
tion is used in this stage. 

The main obstacle to these uses was the size of the available simulations. Although it was to be 
a tool for autopilot research, the simulation program for V/STOL aircraft with rigid-body dynamics 
threatened to be much larger than the autopilot program. However, the effort and labor required 
are proportional to size for a digital program, so it is essential to minimize size, especially in the 
critical early stages of the research. Much was gained by properly partitioning the simulation and 
several partitions that have proved advantageous to our work are described. Since these partitions 
are inevitably motivated by the nature of flight, a description of flight is first given and then the 
system is defined. 


DESCRIPTION OF FLIGHT 


The motion of a rigid aircraft can be thought of in terms of its six degrees of freedom, three of 
which describe the path of the center of gravity through space and three describe the rotational 
motion of the aircraft about its center of gravity. 

The trajectory is defined in some inertial frame. Figure 1 shows a coordinate frame with origin 
at the touchdown point on the runway. This frame is suitable to describe terminal area approach 
paths, and it is an inertial frame if a flat earth is assumed. The trajectory is controlled through the 
aerodynamic and engine forces, which can be described by their body axis components. The body 
axes are a frame attached to the structure with in the aircraft plane of symmetry (fig, 1). If 

the velocity vector is maintained in the plane of symmetry, then the aerodynamic forces as well as 
the engine forces all lie in the plane of symmetry. The aerodynamic force is normally given by its 
components along and perpendicular to the velocity vector (called drag and lift, respectively). This 
force vector is located with respect to body axes by the angle of attack, a, and is a strong function 
of a so that lift force is controlled by controlling a. The other major force is the engine thrust, 
whose magnitude is controlled by a throttle and, for the AWJSRA, its direction with respect to the 
body axes is controlled by a rotatable exhaust nozzle. In addition, compressed air from the engine is 
exhausted through a special flap to augment the aerodynamic lift during landing. Thus, control over 
the inertial path requires control over the inertial orientation of the plane of symmetry as well as 
forces in that plane. 

The inertial orientation is usually defined by the Euler angles <F, 0, and (fig. 2). Control over 
and \p provides control over the inertial orientation of the plane of symmetry and control over 0 
provides control over the angle of attack and therefore the Uft. Control over these attitudes is 
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obtained with the applied moments. For this purpose, small, movable aerodynamic surfaces are 
found at the extremes of the structure (fig. 2) — ailerons at the wing ends and tail and rudder at the 
end of the body — which provide control moments about three orthogonal directions. The terms 8a, 
be, br refer to surface deflection angles. This location at the structural extremes also minimizes the 
disturbance forces associated with attitude control and hence minimizes the coupling of attitude 
controls with the center of gravity motion. 

An experienced passenger may have noticed that there is little rotational motion in commercial 
flights or rather he may have noticed the discomfort when it is otherwise. The attitude is almost 
always steady or slowly varying, and this characterizes the attitude motion generated with the 
simulation. There are also small-amplitude variations about this steady motion, typically at high 
frequency compared to variations in the center-of-gravity motion because the aircraft has inherently 
faster dynamics and controls in the rotational degrees of freedom than in the linear degrees. 


AIRCRAFT SIMULATION - MAIN SUBSYSTEMS 


Figure 3 shows a subdivision of the simulation at its broadest level of subsystems and also the 
input-output relations among them. The autopilot outputs are commands to the six force and 
moment controls of the aircraft model. The aircraft model calculates the applied engine and aero- 
dynamic forces and moments that are output to the equations of motion section which integrates 
these equations and calculates environmental parameters (such as winds, density, etc.). State infor- 
mation and environmental data are then fed back to the autopilot and aircraft model. In the process 
of simulation, the flight evolves by cycling through this loop, where each cycle corresponds to the 
integration step size used in the numerical integration of the equations of motion. 

The principle intent of this modular structure is to minimize the programming labor needed 
for autopilot research by isolating the equations into blocks according to the frequency with which 
they will change. The equations of motion subsystems apply to all aircraft and all autopilots. The 
aircraft model changes rarely as new modeling data are obtained from flight identification, or it 
changes among groups working on different aircraft. Finally, the autopilot changes most fre- 
quently. The labor of plugging in a revised or alternate subsystem is minimized since there are very 
few interface variables. 

Finally, the same structure can be used for all simulations in the anticipated hierarchy of 
simulations so that the autopilot can be transferred up the hierarchy with a minimum of plug-in 
labor, including, eventually, the aircraft flight-control system itself. While the interfaces have been 
minimized, complexity and large size remains internalized in the subsystems, particularly the air- 
craft model and the autopilot. To consider further partitioning of the system and also to define it, 
the equations for each of the three main subsystems are reviewed next. 


Equations of Motion and Environment Subsystems 

Figure 4 shows the equations of motion. Rigid body and flat earth are assumed as is appro- 
priate in our case. These equations exhibit the nature of the dynamics without containing any detail 
about the aircraft or autopilot, so this subset of system equations will appear for any aircraft or 
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autopilot. Where vectors occur in the notation, a subscript indicates the coordinate frame in which 
its components are given, for example, Rj is the position vector in inertial coordinates. 

The first two equations describe the trajectory in the inertial reference frame. The aircraft 
forces are assumed given in body axes and must be transformed to the inertial frame. During a 
simulation run, these are integrated to give the trajectory. The second pair of equations describe the 
rotational degrees of freedom. The angular acceleration equation is written in body axes because the 
aircraft matrix, I, is constant in this frame. The applied moments are assumed given in body axes 
also. 


The system dynamics are partitioned later among the translational and rotational degrees of 
freedom. Attitude enters the trajectory motion through the transformation, and to the extent 
that the force Fg depends on variables of the rotational motion, although this is not shown here. 
The translational degrees enter the rotational motion through only, although this dependence is 
not yet shown. 

Generally, the air mass will be in motion relative to the ground and this enters the problem 
because the aircraft forces are functions of the velocity relative to the air mass rather than the 
inertial velocity. The wind is modeled as the sum of a steady component {WB) and a gust com- 
ponent iWG) generated from a linear system driven by Gaussian distributed noise; this model is 
frequently used in simulations. Parameters of the atmosphere — density, temperature, pressure, and 
speed of sound — are also required for calculating aircraft forces and these are generated from an 
altitude-dependent, standard atmosphere model. 

Finally, various functions of the motion must also be calculated in the simulation: the Eular 
angle transformation and the matrix. A, used in the equations of motion and also the velocity 
vector with respect to the air mass required by the aircraft model. This vector appears in the model 
in the form of its spherical coordinates {VA,oc, |3). Additional quantities used by the aircraft model 
in general are the dynamic pressure, Q and the angle rates a and |3. 

This group of equations is the subset of system equations that define the equations of motion 
subsystem of the simulation. This subsystem appears in any simulation independent of aircraft or 
autopilot. The inputs to this subsystem are the total aircraft forces and moments in body axes and 
the aircraft mass, and its outputs are the translational and rotational motion variables, wind, and 
atmospheric parameters, and various parameters of the velocity vector with respect to the air mass. 
Figure 5 is a block diagram of this subsystem. 


Aircraft Model Subsystems — AWJSRA 

The aircraft model subsystem is defined by its equations (fig. 6). This subsystem is the princi- 
pal source of size and complexity in the simulation and in the control problem so the equations are 
given in sufficient detail to show these aspects. Where details of a specific aircraft occur, these are 
taken from the augmentor wing. The function of this subsystem is to calculate the aircraft forces 
and moments that arise from the engine and from the airflow over the aircraft structure. 

Engine— The engine model requires first that sufficient data be stored to define fuel flow, 
thrust magnitude, and the equivalent augmentation thrust. This last is a fictitious thrust force that 
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measures the amount of compressed engine air exhausted through the trailing edge of the experi- 
mental flap of the augmentor wing. The stored functions are nonlinear functions of atmospheric 
parameters, airspeed, and throttle control. The fuel consumption accounts for slow variations in 
aircraft mass during flight. The engine dynamics that govern variations in thrust in response to 
throttle commands are modeled simply as a first-order system. However, engine dynamics are 
normally quite complex, reflected here in the form of a time constant that changes (depending on 
whether thrust is increasing or decreasing) and in the form of constraints on the extremes of 
magnitude, rate, and acceleration of thrust. The nozzle angle dynamics are taken simply as a 
constant rate. The parameter, Cj, is the dimensionless equivalent thrust and appears as a variable 
throughout the aerodynamic functions. The factor QS is commonly used throughout aerodynamics 
to make forces dimensionless. Finally, the engine force is given by the thrust and directed relative to 
the body axes frame by a unit vector given from the nozzle angle, and the engine moment is given 
straightforwardly (fig. 6). This description is specialized to the AWJSRA, particularly in the appear- 
ance of Cj and v. The term Cj appears throughout the aerodynamics, and its contribution to 
modeling complexity is apparent. The nozzle appears only in the very simple manner shown, but it 
contributes major difficulties to the solution of the control problem, both conceptually and compu- 
tationally. Engine models for other powered lift aircraft would also have special effects that would 
increase the complexity, although different in detail. 

Actuators— The dynamic response of the aerodynamic control surface deflections to control 
commands is considered next (fig. 6). The actual hardware by which the pilot or autopilot com- 
municates with these surfaces is a complex sequence of linkages, servomotors, hydraulic actuators, 
and the surface itself, as defined on an engineering schematic of some kind. However, we are not 
interested in this level of detail, so an input-output model of the actuator system is identified from 
its detailed description and the dominating dynamics separated from that and modeled in the 
simulation as a simple first- or second-order system. The remaining dynamics are much faster than 
those of the aircraft motion. The remaining features of concern are the limits on surface position 
and rate, resulting from aerodynamic or servomotor saturation, which are reflected in the control 
problem as limits on control power. 

Aerodynamics— The principal source of size and complexity in our system is the aerodynamic 
model. Generally, the aerodynamic forces and moments are functions of time histories (up to the 
present) of motion and control variables. However, almost alTflight can be conceived as being small 
departures from some steady flight condition and this motivates modeling the aerodynamic force as 
a sum of terms (fig. 6). The leading term,/^, is the force in steady flight and the remaining terms 
are superposed independent small effects that account for nonzero values of various motion and 
control variables. For completeness, the expression is premultiplied by a transformation 
from wind axes components in which the aerodynamic force data are usually given, and the factor, 
QS, that separates the principal dependence of forces on air speed and also makes the vectors in the 
sum dimensionless. 

The leading term,/^, is a function of a, F, Cj, and altitude, and it accounts for the principal 
dependence of aerodynamic forces on these variables. Although the terms are functions of 
several variables, the principal variations in flight and their difference for zero results from the 
variable 5z, which is one of the angles ||3, 5a, be, 5i?[. These terms are either linear or nearly so in 
the variable bi. For the controls, their usefulness is closely dependent on their near-linearity. The 
remaining terms F,1 account for various angular rates ja, |3, whose effects are modeled as 
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linear and are made dimensionless by the appropriate factor from {My*} • Similar remarks apply to 
the description of the moments. 

The components of the force and moment vectors given in wind axes are individually identi- 
fied in the working nomenclature as coefficients of drag, Y, and lift forces and as coefficients of 
roll, pitch, and yaw moments — one coefficient for each component of each vector in the sum. 
These coefficients are functions of varying degrees of complexity and sufficient data to define all 
the coefficients in the list are stored in the simulation. The gross properties of this collection of 
coefficients are examined next. 

Aerodynamic coefficients— Figure 7 summarizes in a matrix the coefficients required for the 
aerodynamic mode. The columns itemize the terms that are summed to form the total force and 
moment and the rows itemize the wind axes components of each term. The resulting nomenclature 
of the coefficients thus identifies each contribution to aerodynamic forces and moments by source 
and by effect — this is the working nomenclature. A large body of classical literature exists which 
analyzes the origins and consequences to the flight dynamics of each coefficient. 

The collection of coefficients divides about evenly into trivial and nontrivial ones. The non- 
trivial ones vary in complexity from constants to functions of four variables. Much of this com- 
plexity is characteristic of powered lift aircraft compared to conventional aircraft, for example, the 
engine parameter, Cj, appears in most of the coefficients. And parameters such as F, 5a, which 
conventionally have linear additive effects, are now nonlinear and even nonadditive as for F. 

In the simulation program, these coefficients are evaluated at each integration step and this 
constitutes a large part of the computational and storage load in a six-degree-of-freedom simulation. 
This situation is common to powered lift vehicles, and efficient computational methodology is a 
significant programming problem that is also reflected as a significant programming problem in the 
autopilot. The functions are handled by storing tabulated functions that are then interpolated 
during a simulation run. Figure 8 shows an example — a lift-drag polar that defines the leading force 
term by mapping values of the three independent variables (a, F, Cj) into values of Cj^ and Cyy. For 
a conventional aircraft, we would have a single independent variable and therefore a single line on 
one plot. For the AWJSRA, two additional variables require first additional lines and then addi- 
tional plots. 

Half of the 60 possible coefficients are zero. This reflects the physical symmetry of the aircraft 
as well as the normal operation of the aircraft with the velocity vector in the plane of symmetry. If 
the rows and columns of the matrix of the coefficients are shuffled, the structure shown in figure 9 
appears. As seen, the system of variables and of force and moment components almost divides into 
two independent subsystems (referred to as the longitudinal and lateral degrees of freedom), which 
separates motion in the plane of symmetry from motion lateral to it. If we were to linearize the 
equations of motion about an unaccelerated flight condition, this decomposition appears through- 
out the systems of equations, and this has been exploited heavily in the classical development of 
aircraft control since the small-disturbance behavior of the aircraft can thus be analyzed as two 
subsystems half the size. 

This subdividion is now exploited in the nonlinear simulations, nor in our control problem, 
where we are interested in the inertial path of the vehicle as much as its local small-perturbation 
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dynamics. It wDI prove more rewarding in terms of solving the autopilot problem to exploit a 
division into linear and rotational degrees of freedom. Nevertheless, this shows that there are other 
partitions of the complete system into loosely coupled subsystems besides the ones to be focussed 
on here. 

Figures 6 through 9 have defined the model of the aircraft forces and moments that result 
from the engine and aerodynamics. Their interrelationships are exhibited in a block diagram of the 
aircraft model (fig. 10). The internal blocks correspond to the main subroutines of the digital 
simulation of this subsystem; the engine is modeled in a single block while the aerodynamic 
description is divided into three blocks — actuators, forces, and moments. 

The input-output relations among these blocks are arranged as two main parallel streams of 
cause and effect. This divides the calculations into quantities for which significant variations occur 
slowly or rapidly in time relative to each other. The lower stream contains slowly varying quantities 
— engine model, aerodynamic forces, and many coefficients of the moment equations which vary 
slowly with a, F, and C.-. Most calculations of the aircraft model can be separated into this stream. 
This partitioning is exploited when computation time is critical by carrying out the slow stream at, 
for example, half the frequency of the fast stream. This has been done in simulations operated in 
real time, primarily man-in-the-loop simulations. On the one hand, the computation time fora cycle 
must not exceed the corresponding real time increment being simulated, and that time increment 
should also be about 10 times the natural frequency of the dynamics being simulated. This parti- 
tioning of calculations was used when the computational load required for simulating powered lift 
vehicles exceeded the computer speed if carried out in a single stream. An alternate solution, also 
used, is a faster computer. A similar constraint on computational cycle time will occur for the 
autopilot in flight, and a similar partitioning of the control calculations into subsystems evaluated at 
different rates can be exploited there as well. 

Autopilot— In the simulation of aircraft flight, the motion is always controlled, that is, the 
values of all controls must somehow be defined to operate the simulation. In the present case, 
where the simulation is used for autopilot development, experimental autopilot logic fills this role. 
Although the autopilot is not discussed in detail, it is useful to describe its functioning at a general 
level since this will motivate a useful partitioning of the aircraft simulation (fig. 1 1 ). 

The purpose of the autopilot, in general, is to fly some reference path defined, for example, by 
air traffic control. To do this, the autopilot necessarily solves the equations of motion for the 
control settings that fly the reference path and corrects any errors from the path. 

Figure 1 1 shows the input information to be the state variables divided into variables of the 
translational and rotational degrees of freedom, and also wind information and inputs from air 
traffic control that define the reference trajectory. The outputs are commands to the control 
actuators, subdivided into force and moment controls. 

Internally, the autopilot logic is partitioned into trajectory control and attitude control. The 
attitude control logic calcualtes the moment control commands that execute attitude and angular 
velocity commands from the trajectory control logic. The trajectory control logic uses the force 
controls, UF, to control forces in the plane of symmetry; 5th and v are engine controls and F, the 
flap deflection angle, controls lift through its dependence on wing configuration. It also uses the 


172 



attitude control subsystem to control the orientation of the plane of symmetry in inertial space and 
it uses pitch angle to control aerodynamic force through angle of attack. 

The control problem is thus divided into two smaller problems based on a partitioning of the 
system into translational and rotational degrees of freedom. The trajectory is controlled by large 
engine and aerodynamic forces that inherently have slower dynamics than the control moments 
used by the attitude control subsystem. In addition, it is assumed that attitude commands from the 
trajectory control logic will have only low-frequency variations for the passenger-oriented opera- 
tions considered here. Thus, for the controlled airplane, the rotational subsystem appears to the 
translational subsystem as a fast servo insofar as attitude response to attitude commands is con- 
cerned. This result can be used to simplify the rotational degrees of freedom and to obtain a 
relatively small simulation that has proved extremely useful in studying the trajectory control 
problem (fig. 1 2). 


SIMPLIFIED SIX-DEGREE-OF-FREEDOM SIMULATION 


Figure 12 shows the complete system — autopilot, aircraft model, and equations of motion — 
subdivided into translational and rotational subsystems that separate the trajectory variables and 
attitude variables. The translational subsystem is affected by the rotational subsystem through the 
attitude response to attitude commands and through a perturbation force that results from high 
frequency rotational dynamics and has negligible effect on the trajectory. Because of the faster 
rotational dynamics, the complete rotational subsystem can be removed and its input-output 
behavior obtained approximately from a linear second-order system (fig. 1 2). 

Thus the simulation can be partitioned into a half-size simulation to study the trajectory 
control problem. This has proved an ideal tool, much smaller in size and computation time than the 
parent six-degree-of-freedom simulation. In fact, the size of the trajectory control logic program is 
twice that of the remainder of the simulation, which gives a comfortable ratio of the tool size to 
problem size for this portion of the research. If the complete six-degree-of-freedom simulation were 
to be used for the trajectory control research, then the reverse would be true, and the control logic 
program would be outweighed by the rest of the simulation by more than two to one. 


HIERARCHY OF SIMULATIONS 


Figure 13 shows some aspects of literal size revealed by the hierarchy of simulations. A 
six-degree-of-freedom simulation was available for our research needs, but the size of this object was 
a problem, as well as its organization, which was oriented toward general-purpose use with man-in- 
the-loop simulation. It was large in terms of the number of boxes of FORTRAN cards that could be 
carried about lightly, or the number of listings that could be spread out on a desk, or the number of 
sheets of paper required to write down a working definition of the system. In this regard, the report 
that defines the aircraft model has 96 pages of equations and tabulated function data (ref. 3). These 
sizes are significant in relation to the desire for easy accessibility, use, understanding, and 
modification during the critical early stages of the research when the central conceptual problems 
are solved. 
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In response to this problem of size, a hierarchy of simulations of increasing size and detail has 
been assembled or planned. First, the simplified six-degree-of-freedom simulation described in 
figure 1 2 is used to develop the trajectory control logic. The shaded and unshaded portions of the 
bar in figure 1 3 are the size of the trajectory control logic and the size of the remainder of the 
simulation, respectively, in feet of FORTRAN cards. Not only is the absolute size of this program 
quite small, but the size of the research portion of the simulation is twice that of the rest of the 
simulation program. This is a comfortable working ratio between the size of the tool used and the 
size of the problem being solved. Although the simulation is very much smaller than the six degrees 
of freedom, the control problem remains difficult. 

The second member of the hierarchy is the six-degree-of-freedom simulation, used for attitude 
control and navigation studies, as well as for verification of the trajectory control logic in a more 
detailed simulation of the flight dynamics. This tool is much larger and the ratio of tool size to 
problem size is less favorable. (One response to tool size is to enlist a specialist to work the tool for 
you, but at the expense of some loss of familiarity with the system being controlled.) 

The third member of the hierarchy of simulations is a large man-in-the-loop laboratory simula- 
tion on computers dedicated to such simulations. The added detail included a cockpit, displays, 
sensors, and also a copy of the STOLAND airborne computer and peripheral equipment that may 
be used for a thorough evaluation before any flight test. With a simulation of this size, a staff of 
specialists and assistants who work solely with the simulation is required. 

The organizational size itself creates certain problems of inconvenience to the research. The 
demand for the simulator-STOLAND is sufficiently great that long lead time scheduling or incon- 
venient working hours become necessary and it is quite possible that an additional simulation, 
intermediate in size between the second and third members of the hierarchy, will be programmed to 
minimize the time needed on simulator-STOLAND. 


SUMMARY 


Simulations of V/STOL aircraft for autopilot research were introduced as an example of a 
large-scale system. Large, in this case, was perceived in relation to the uses of the simulation, that is, 
in terms of the research requirement for a simulation that was easy to work with and in terms of the 
complexity and difficulty of the control problem. Both problems were attacked and at least par- 
tially solved through partitioning of the system of equations and hierarchic organization of the 
system. The various organizations described are listed in figure 14. 

First, a hierarchy of simulations was assembled — a sequence of increasingly detailed simula- 
tions of the aircraft paralleling a corresponding hierarchic development of the autopilot. The object 
was to have a digital simulation no larger than necessary at each stage of the research since such 
factors as ease of access, use, understanding, and modification are all proportional to size. 

Second, a modular organization was given to the simulations, which were divided into three 
main subsystems — the autopilot, the aircraft model, and the equations of motion and environment 
section. These modules have a minimum of interface variables and provide for plug-in simplicity so 
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that frequent revisions of the autopilot and occasional changes of aircraft occurring in the research 
can be made. 

Third, the dynamics of the system were subdivided into translational and rotational degrees of 
freedom, based on the different frequencies at which significant variations in motion variables and 
control forces and moments occur in the two subsystems. This is the basis of the autopilot parti- 
tioning into two smaller control problems. This frequency separation can also be used to satisfy the 
real-time constraint on the autopilot cycle time should the computational load for powered-lift 
aircraft prove too large in relation to computer speed. This was also the basis for devising a 
simplified six-degrce-of-freedom simulation by replacing the rotational subsystem with its input- 
output relation. This has been a very convenient research tool for the trajectory control problem 
and has permitted a study of that problem independent of the attitude control problem. 

Fourth, an alternate partitioning of the aircraft dynamics into longitudinal and lateral degrees 
of freedom was mentioned. This is used in the classical linearized analysis and separates the 
dynamics into motion in the plane of symmetry and lateral to it. 

Finally, partitioning the aircraft model section into two subsystems whose variations were well 
separated in frequency permitted these subsystems to be evaluated at different rates during simula- 
tion. This has been used in real-time simulations to meet the real-time constraints on computation 
cycle time, especially for powered lift V/STOL aircraft. 
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Figure 1.— Translational degrees of freedom. 
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Figure 2.— Rotational degrees of freedom. 
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Figure 3.— Block diagram of AWJSRA aircraft simulation. 



Figure 4.— Equations of motion — rigid body, flat earth. 
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Figure 5.— Block 



R|,v,,v, 

ATT, fig 
V ,13 ,ot ,0 
W|,AT 


of equations of motion and environment model. 
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Figure 9 — Aircraft model — coefficients for longitudinal and lateral degrees of freedom. 


GROUT: (R|, V,, fig, VA, a, P, a, p, AT, . . .) 



Figure 10.— Block diagram of aircraft model — (AWJSRA). 
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Figure 12.— Subsystems — translational and rotational degrees of freedom. 
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Figure 13.— Hierarchy of simulations. 


1. 

HIERARCHY OF SIMULATIONS 

LIMIT SIZE OF SIMULATIONS 

2. 

MODULAR ORGANIZATION 

plug-in SIMPLICITY FOR AUTOPILOT 
OR AIRCRAFT CHANGES 

3. 

TRANSLATIONAL/ ROT AT I ONAL 
DECOMPOSITION 

SUBDIVIDE CONTROL PROBLEM 
REDUCE COMPUTATION CYCLE TIME 
SIMPLIFIED SIX-DEGREE-OF-FREEOOn 
SIMULATION 

4. 

LONGITUDINAL/LATERAL 

DECOMPOSITION 

SUBDIVIDE LINEARIZED DYNAMICS 

5. 

AIRCRAFT MODEL "SLOW'V'FAST" 
CALCULATIONS 

REDUCE COMPUTATION CYCLE TIME 


Figure 14.— Summary of organizations. 




